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High resolution measurements of superfluid density ps{T) and broadband quasiparticle conduc- 
tivity cri(Q) have been used to probe the low energy excitation spectrum of nodal quasiparticles in 
underdoped YBa2Cu306+y. Penetration depth A(T) is measured to temperatures as low as 0.05 K. 
(TiiQ) is measured from 0.1 to 20 GHz and is a direct probe of zero-energy quasiparticles. The 
data are compared with predictions for a number of theoretical scenarios that compete with or 
otherwise modify pure dx2_y2 superconductivity, in particular commensurate and incommensurate 
spin and charge density waves; d^2_y2 + is and d^2_y2 -\-idxy superconductivity; circulating current 
phases; and the BCS-BEC crossover. We conclude that the data are consistent with a pure d^2_y2 
state in the presence of a small amount of strong scattering disorder, and are able to rule out most 
candidate competing states either completely, or to a level set by the energy scale of the disorder, 
~ 4 K. Commensurate spin and charge density orders, however, are not expected to alter the 
nodal spectrum and therefore cannot be excluded. 



I. INTRODUCTION 

The physics of the cuprate high temperature super- 
conductors is that of strong Coulomb repulsion in nearly 
half- filled Cu02 planes As charge carriers are doped 
into these materials, the two most prominent electronic 
states are the antiferromagnetic (AFM) Mott insulator 
and the d-waye superconductor. While the AFM and 
the optimal-to-overdoped superconductor appear to be 
well understood, the physics of the underdoped part of 
the phase diagram that lies between them remains firmly 
incompatible with standard theory. The most prominent 
feature of this region is a pseudogap that suppresses low 
energy spin and charge fluctuations and persists above 
the superconducting transition to a temperature T* 
The pseudogap temperature is highest close to the Mott 
insulator and decreases monotonically as doping, p, is in- 
creased towards optimal doping. Identifying the nature 
of the pseudogap state remains a difficult and open prob- 
lem. 

States of matter are characterized by their symme- 
tries and their low energy excitation spectra, d-wave su- 
perconductivity, for instance, breaks four-fold rotational 
symmetry and is distinguished by the presence of nodal 
quasiparticles with a characteristic linear energy spec- 
trum. The (i-wave state in the cuprates was first iden- 
tified from observations of a linear temperature depen- 
dence of penetration depth. A, and superfluid density, 
ps = 1/A^.^ The ability of superfluid density to couple 
directly to itinerant electronic degrees of freedom gives it 
the potential to be a sensitive thermodynamic probe of 
pseudogap physics, with many candidate states expected 
to leave characteristic signatures in the low energy quasi- 
particle spectrum. Here we search for these signatures 
using high resolution measurements of penetration depth 
and broadband quasiparticle conductivity, made on very 
clean crystals of underdoped YBa2Cu306+2y 



There have been a wide range of proposals put for- 
ward to explain the cuprate pseudogap. In one impor- 
tant category, strong pair correlations are already built 
into the normal state. This scenario has its roots in 
Anderson's resonating- valence-bond spin liquid,^ and the 
idea that pair correlations emerge directly from the Mott 
insulator remains a compelling proposition. The 'gos- 
samer superconductor' — a BCS wavefunction in which 
double occupancy has been heavily suppressed — typ- 
ifies this approach and may provide a useful represen- 
tation of the underdoped electronic state. ^ The impli- 
cation for the phase diagram is that T* marks the for- 
mation of tightly bound Cooper pairs, with low phase 
stiffness and strong quantum and thermal phase fluctu- 
ations heavily suppressing Tc^^^^^^^^ At tempera- 
tures not too far above the superconducting transition, 
the idea of pre-existing pairs finds support from a num- 
ber of experiments: terahertz spectroscopy reveals a fi- 
nite phase-stiffness;^^ Nernst-effect measurements appear 
to detect the phase-slip voltage from thermally diffus- 
ing vortices fiSiii high-field magnetometry reveals excess 
diamagnetism^ and STM^^ and /iSR32. detect what ap- 
pear to be droplets of precursor superconductor. Re- 
lated to this, the theory of the BCS to Bose-Einstein 
condensate (BEC) crossover makes a prediction that can 
be tested here: a T^/^ power law in ps{T)^ due to the 
direct thermal excitation of bound Cooper pairsJ ^'^^ 

Another class of proposals seeks to explain the pseudo- 
gap in terms of competing orders and quantum criticality. 
In such a scenario, T*(p) marks the boundary of a dis- 
tinct thermodynamic phase; must be accompanied by a 
broken symmetry; and goes to zero at a quantum crit- 
ical point within the superconducting phase. This idea 
was initially motivated by the observation near optimal 
doping of so-called marginal Fermi liquid behaviour^ 
in which unusual power laws in resistivity p(T), optical 
conductivity (Ji{Q) and other physical quantities could be 
understood in terms of scattering from a scale- invariant 
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fluctuation spectrum, as would be expected near a zero- 
temperature critical point On crossing T*, these fluc- 
tuations should generically condense to form the broken 
symmetry state of the pseudogap phase. While there is 
evidence of an AFM quantum critical point in electron- 
doped materials, the situation on the hole-doped side 
is much less clear. Identiflcation of a particular compet- 
ing order that appears at T*(p) would have strong im- 
plications not just for the pseudogap, but for the origin 
of non- Fermi- liquid behaviour elsewhere in the cuprate 
phase diagram. 

Competing orders are in fact prevalent in the cuprates, 
in part as a result of the extreme sensitivity of the 
doped Mott insulator to perturbations^ Outside the 
AFM phase, long-range magnetic order is replaced by 
glassy spin correlations ^2L2^i2^ although this short-range 
magnetism is likely a response to chemical disorder ^^i^ 
Neutron scattering experiments on La2-a:Sra^Cu04 have 
revealed incommensurate spin correlations in supercon- 
ducting samples^i^i^ that were later identifled as stripe 
ordering of spins and holes.— ^ Stripe correlations ap- 
pear to be widespread in the underdoped cuprates, and 
are particularly strong near p = | doping. In ap- 
plied fleld, the suppression of superconductivity in vortex 
core^i^i leads to co-existing superconductivity and spin- 
density- wave order i^i^i^i^i^ Scanning tunneling spec- 
troscopy of the vortex cores reveals that this is accompa- 
nied by prominent checkerboard charge-density order 4^ 
Similar four-lattice-constant modulations of the density 
of states are seen in zero fleld at various points in the 
phase diagrami^i^ For the most part, these competing 
orders occur in narrow ranges of doping; or in particular 
materials; or in response to external perturbations such 
as point disorder or applied magnetic fleld. While they 
attest to the complexity of the doped Mott insulator, 
they offer only hints at the physics underlying the forma- 
tion of the pseudogap. 

The lack of compatibility of the observed ordered states 
with T*(p) has led to interest in 'hidden orders', in 
which the broken symmetry is subtle and difficult to 
detect with standard scattering experiments. Proposals 
include circulating current phases that preserve transla- 
tional symmetry^i^i^ and orbital antiferromagnetism, 
for example the (i-density wave state (DDW)."^^ Inter- 
estingly, a set of recent experiments now appears to 
have detected signatures of one or more of these phases. 
^SR^ and polar Kerr effect^^ have established the on- 
set of time-reversal-symmetry breaking (TRSB) at T*(p) 
in YBa2Cu3 06+2/7 but the signals are extremely weak. 
It has been suggested that a variant of the DDW, the 
dxy + idj;2_y2 density wave^^ would contain a subdom- 
inant but macroscopic TRSB component and be consis- 
tent with the small magnitude of the observed effects. 
Spin-polarized neutron scattering on YBa2Cu306+2y has 
detected weak signatures of a novel magnetic order that 
preserves translational symmetry^ and has a form con- 
sistent with the 0/7 circulating current phase proposed 
by Varma^ shown in schematic form in the inset of 



FigO The detailed picture is complicated by the pres- 
ence of an in-plane component of magnetic moment, al- 
though it has been suggested that this could arise from 
orbital currents that circulate through apical oxygens 
while preserving the B// symmetry. It also remains to 
be seen how ubiquitous the effects are: /iSR experiments 
on La2-ccSra^Cu04 have so far failed to observe TRSB^ 
but have not yet been carried out with the same sen- 
sitivity as Ref. JEdu In contrast, new neutron scattering 
experiments^^ on IIgBa2Cu04+(5 have detected the same 
type of 0/7 magnetic order seen in YBa2Cu306+2y 
we will discuss in more detail below, this type of order 
has a strong effect on the low energy states of the super- 
conductor, and should be highly visible in measurements 
of Ps{T). 

Finally, there have been suggestions that pure (1^2 _y2 
superconductivity may compete with superconducting 
states of different symmetry^i^i^ motivated in part by 
reports of anomalously large inelastic scattering of nodal 
quasiparticles below Tc^^^ This critical-like scattering 
has been shown to be compatible with a quantum phase 
transition to a dx2_y2-\-is or dx2-y2-\-idxy state. To date, 
there is a limited amount of direct experimental evidence 
in support of such phases^i^i^ — here we use measure- 
ments of superfluid density to place tight constraints on 
the existence of such states in YBa2Cu306+2y 

This paper is organized as follows. In Sec. [Til we show 
how measurements of penetration depth and broadband 
microwave conductivity can together be used as a probe 
of the quasiparticle excitation spectrum and the struc- 
ture of the superconducting energy gap. In Sec. [TTIl we 
catalog how different competing orders affect the super- 
fluid density, including the effect of disorder. In Sec. IIVI 
we introduce the experimental methods used to measure 
superfluid density and broadband microwave conductiv- 
ity. Results are presented and discussed in Sec. fol- 
lowed by a summary of our conclusions in Sec. IVI[ Ap- 
pendix [A] presents analytic results for the effect of disor- 
der on dx2_y2^ dx2_y2 -\-idxy and dj;2_y2 + is-type super- 
conductors with isotropic Fermi surfaces, and shows how 
this eventually blurs the distinction between dx'^_y2 and 
dx2-y2 -\- idxy states. 



II. PENETRATION DEPTH AND MICROWAVE 
CONDUCTIVITY 

Microwave experiments can be used to probe the low- 
energy excitation spectrum of a superconductor in two 
ways: through the temperature dependence of the pene- 
tration depth A(T); and from broadband measurements 
of the oscillator strength in the flnite-frequency quasi- 
particle conductivity spectrum (Ji{Q^T). The theory of 
penetration depth and microwave conductivity of uncon- 
ventional superconductors has been developed in great 
detail^^^^^^^^^ but useful insights about low- 
lying excitations can be obtained from the weak-coupling 
BCS theory. For the case of an isotropic Fermi surface. 
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which should be adequate for describing the low-lying 
excitations in the cuprates, A(T) is given b y^^i^^ 

= ^ du tanh — — ^ . (2) 

Here Aq is the zero-temperature penetration depth 
in the absence of disorder and competing phases, 
and f((jj/T) is the Fermi function. A Sommerfeld 
expansion reveals the direct connection between A(T) 
and the normalized density of quasiparticle states 
N{uo): if N{uo) = A^o + Niuo + \N2Uo'^ + ... then 
\l/\\T) = 1 - TVo - 21n27Vi^BT - -^N2{kBTf - ... 
The residual density of states (DOS) Nq represents 
zero-energy excitations, which arise in a superconductor 
either from impurity pair-breaking, or from certain types 
of competing order, notably the 6//-type circulating 
current phase^*^. Note that Nq does not appear in the 
temperature dependence of A, but instead results in a 
deviation of A(T^O) from Aq. This shift in penetration 
depth is difficult to resolve experimentally, because 
Ao is neither known a priori^ nor can the absolute 
value of A(T 0) usually be measured with sufficient 
accuracy. However, a direct determination of Nq can be 
obtained from the uncondensed weight spectral in the 
quasiparticle conductivity (Ji(l^,T). From the oscillator 
strength sum rule, 

No = l^^\l I ai(l],T^O)dl], (3) 
Jo 

where Q^c is a frequency cut-off chosen to capture the 
oscillator strength of the conduction electrons only. 

In a clean-limit BCS superconductor, N{(jo) is deter- 
mined by the k-space structure of the superconduct- 
ing order parameter Ak: N{ijo) = Re(a;/ ^/o;^ — A^)pg, 
where (•.•)pg denotes a Fermi surface average. This 
makes psiT) a sensitive probe of order parameter sym- 
mery. In particular, for a d-wave superconductor in 
two dimensions, the linear dispersion of Ak about the 
gap nodes leads to N{(jo) cx uo and Aps{T) ex T. An 
s-wave superconductor, by constrast, usually has a fi- 
nite energy gap and shows activated behaviour, ps{T) oc 
exp(— Amin/^B^), where Amin is the minimum of the 
energy gap on the Fermi surface. The effect of impu- 
rity scattering on N{u;) and psiT) is important and is 
reviewed in Appendix [Aj where we give analytic results 
for dx2_y2^ dx2_y2 -\-idxy and dx2_y2 -\-is superconductors 
with isotropic Fermi surfaces in the presence of point de- 
fects. The main effect of disorder is for the quasiparticles 
to acquire a lifetime, the magnitude and energy depen- 
dence of which depend on the concentration and scat- 
tering strength of the defects. Near the unitarity limit, 
scattering leads to a zero-energy resonance that overlaps 
with the continuum of quasiparticle states in the dj.2_y2- 
wave superconductor, resulting in a residual density of 



states in N{uj) and a crossover to behaviour in ps{T). 
This also happens for the dx2_y2 + idxy superconductor, 
despite there initially being a finite gap in the excitation 
spectrum. As a result, above a certain level of disorder, 
dx'2-y2 and dj;2-y2 -\-idxy states become impossible to tell 
apart using microwave spectroscopy. In Fig. [9] we show 
how the distinction is lost when the energy scale of the 
disorder, ksTd > A^^^^. The dx'2-y2 + is superconductor 
is different in this respect: nonmagnetic scatterers do not 
cause pair breaking at low energies, and the gap in the 
spectrum is robust. 



III. SUPERFLUID DENSITY AND 
COMPETING ORDERS 

As a sensitive thermodynamic probe that couples di- 
rectly to current-carrying excitations, measurements of 
superfluid density are well suited to detecting changes in 
the nodal quasiparticle spectrum arising from competing 
orders and other physics. A number of authors have in- 
vestigated these effects theoretically. Sharapov and Car- 
botte have performed calculations for a dx2-y2 + idxy 
order parameter and for incommensurate spin density 
waves that nest the nodal points (nested SDW), obtain- 
ing analytic results for ps(T 0) and its leading temper- 
ature corrections.^^ In the absence of disorder they find 
that both the nested SDW and the dx'2-y2 + idxy super- 
conductor have a finite gap everywhere on the Fermi sur- 
face, leading to activated exponential behaviour psiT) ~ 
exp(— Ay/^B^), where A' is the magnitude of the SDW 
or dxy gap. However, nested SDW orders compete for 
Fermi surface, removing nodal states from the T = con- 
densate. In contrast, a transition to a clean dx2-y2 -\-idxy 
state leaves ps{T 0) unchanged. Unfortunately, this 
distinction is difficult to detect experimentally, for rea- 
sons discussed in Sec. [Ill In the presence of disorder, 
both the nested SDW and dx^-y^ + idxy states develop 
a leading quadratic temperature dependence, ps ~ T^, 
similar to that of a dirty dx2-y2 superconductor. How- 
ever, an experimentally detectable difference now arises: 
pair-breaking in the dx2-y2 -^idxy state is accompanied by 
zero-energy quasiparticles, whereas the disordered SDW 
continues to remove low energy states without creating 
a residual DOS. Atkinson has studied the competition 
between nested, incommensurate SDW and dx2-y2 su- 
perconductivity numerically and finds broadly similar 
results^ pointing out that on the basis of the temper- 
ature dependence of ps alone, the effect of disordered 
magnetism cannot be distinguished from dirty but pure 
d-wave superconductivity. He shows that the suppres- 
sion of zero-temperature superfluid density in the nested 
SDW case arises because nodal Cooper pairs cease to 
carry a well-defined current. Modre et al. have stud- 
ied the dx'2-y2 + is pairing state, which also has a fi- 
nite energy gap and activated behaviour in ps{T) at low 
temperature^^ In contrast to the dx'2-y2 -\- idxy case, the 
dx^-y^ -\-is gap is stable in the presence of disorder of any 




FIG. 1: (color online). The presence of a perturbation of the 
form Eq. 2] from the B/j-type circulating currents shown in 
Fig. [21 modifies the nodal spectrum of the (1^2 _y2 supercon- 
ductor in a characteristic way: one node is shifted up in energy 
by 4Acc, one is shifted down, and two are unperturbed. 
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FIG. 2: Individual nodal contributions to the density of states 
N{ijj) from a circulating current perturbation of the form 
Eq. m Inset, upper left: the G/j circulating current pattern 
proposed in Ref. 48. Inset, lower right: an equivalent current 
pattern within a one-band model of the Cu02 planes.— 



strength. In Appendix [A] we show how this arises from 
impurity renormalization of the s-wave gap component. 

Berg et al. have studied the stability of the nodal quasi- 
particle spectrum in the presence of commensurate com- 
peting orders of all types For commensurate perturba- 
tions that do not nest the nodal points, they prove that if 
the perturbation is invariant under time reversal or time 
reversal followed by a lattice translation, the nodal spec- 
trum is stable. While it remains uncertain whether the 
converse holds in general, they examine several important 
cases in which the nodal spectrum breaks down, includ- 
ing certain stripe-like arrangements of spin and charge 
density, and the 9// circulating-current phase that has 
been detected by neutron scattering in YBa2Cu306+2y 
and HgBa2Cu04+5. Confining themselves to a one-band 
model of the Cu02 planes. Berg et al. have used the 
simpler arrangement of orbital currents shown in Fig. [21 
which is equivalent to the 6// state in the more compli- 
cated three-band Cu-0 lattice of Ref. 48. For a pertur- 
bation to the pure d^-^^y^ superconductor of the form 




where rirr' = ±1 is determined by the direction of the 
bond currents in Fig. [21 they find excitation energies 
= + 2Acc {sin(A:a;a) - s\ii{kya) + s\n[{ky - kx)a]}. 
Here is the unperturbed d-wave spectrum and a is 
the lattice spacing. The perturbed nodal spectrum for 
the B// state is plotted in Fig. [H The effect of the 
circulating currents is similar to the Doppler shift from 
a uniform current applied along a diagonal direction: 
one node shifts up in energy by ~ 4Acc, one node shifts 
down, and two are unperturbed. The individual and 
combined contributions to the low energy DOS are 



plotted in Fig. [21 The net effect on N{(jo) is a finite 
residual DOS ^ 2Acc/Ao, and a kink dX uo ^ 4Acc above 
which the linear energy dependence doubles in slope. In 
the clean limit, the superfluid density can be obtained 
from Eqs. [Hand [21 and is plotted in Table [D The limiting 
low temperature behaviour of ps{T) is linear, arising 
from excitations near the two unperturbed nodes. At a 
temperature of order 4Acc//^b, Ps{T) crosses over to a 
second linear regime in which all four nodes contribute 
and the temperature slope doubles. In a clean sample, 
the combination of a residual DOS and a kink in psiT) 
separating two linear regimes should be easily observable 
in experiments. Calculations in the presence of disorder 
have not been carried out, but we expect strong scatter- 
ing impurities to induce additional residual DOS and to 
cause a crossover to behaviour in ps{T)^ as is seen 
in dx2_y2 and dx2_y2 + idxy superconductors. Although 
disorder will mask the effect of circulating currents when 
the crossover temperature > 4Acc/^b, it is expected 
that tight limits on the size of Acc can nevertheless be 
placed, either using Ps{T) or from the magnitude of the 
uncondensed spectral weight in (Ji{Q^ T ^ 0). 

The effect of competing orders on pg (T) and the resid- 
ual DOS is summarized in Table [H The SDW results are 
for the case of ordering wavevectors that nest the nodal 
points. The response to nested charge density waves is 
expected to be broadly similar, with the opening of a 
finite nodal gap that competes for Fermi surface. 



IV. EXPERIMENT 

Measurements of ps{T) and (Ji{Q^T) have been made 
on a single-crystal ellipsoid of YBa2Cu306.333, prepared 
as described in Ref. 7g. Following high pressure anneal- 
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d^2_y2 G/i current loops (1^2 _y2 + is (1^2 _y2 -\-\dxy nested SDW 




777 7/7" 7/7, 7/7, 7/7, 

c c c c c 



residual clean No Yes No No No 

DOS dirty Yes Yes No Yes No 

lowT clean T T/2 ^-A./fc^T ^-A^^^/fc^T ^^-AsDw/fc^T 

Ap.(T) dirty e"^^/'^^^ 



TABLE I: (color online). Effect of competing orders on the superfluid density of a d^2_y2 superconductor. The first column 
shows results for the pure d^2_y2 state. Subsequent columns show the effect of competition from: B//-type circulating 
currents j^^i^^ d^2_y2 + is superconductivity*^ dx2_y2 + idxy superconductivity;^^ and spin density waves (SDW) that nest 
the nodal points.— The first row shows clean-limit excitation spectra for the near-nodal quasiparticles. The second row gives 
the density of states (DOS) N((jj), both for clean systems and in the presence of disorder. Note that the effect of disorder has 
not been calculated for the B/j-type perturbation, and that nonmagnetic disorder has essentially no effect on the d^2_y2 + is 
superconductor. The third row shows the temperature dependence of the superfluid density ps(T), including deviations from 
full condensation as T ^ due to the presence of zero-energy quasiparticles. The fourth row indicates whether a residual 
density of states is expected to be seen in cri(Q, T ^ 0), and the fifth row gives the leading low temperature behaviour of the 
superfluid density. Details of the calculations are given in Appendix [A] for the dx2_y2, d^2_y2 + is and dx2_y2 + idxy states. 
Dirty limit calculations have been made for unitarity limit scatterers. 



ing under a hydrostatic pressure of 35 kbar, controlled 
relaxation of oxygen order in the CuO chains has been 
used to continuously tune Tc in the range 17 K to 3 K. 
Broadband microwave spectroscopy was carried out early 
in the sequence, for Tc = 15.6 K. Measurements of Ps{T) 
in the milliKelvin range were made in the fully relaxed 
state, where Tc = 3 K. 

Ps{T) is obtained from 2.64 GHz surface impedance 
measurements, as described in Refs. [tqI and [80|- The 
sample is positioned at the i^-field antinode of the TE015 
mode of a r utile dielectric resonator, with the microwave 
i^-field oriented along the c axis of the ellipsoid to in- 



duce a6-plane screening currents. Surface impedance 
Zg = Rs -\- iXs is obtained using the cavity perturbation 
approximation: 

+ iAXs = r {A/b(T) - 2iA/o(T)} , (5) 

where A/^(T) is the change in bandwidth of the TEoi(5 
mode upon inserting the sample into the cavity; A/o(T) 
is the shift in resonant frequency upon warming the sam- 
ple from base temperature to T; and F is an empirically 
determined scale factor. The absolute reactance is set 
by shifting AXg (T) so that it matches Rs (T) in the nor- 
mal state. We expect local electrodynamics to be a good 
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FIG. 3: (color online). a6-plane superfluid density ps{T) = 
A^(T) shown at 13 of 37 dopings measured in this study. The 
straight lines are linear fits the data between 5 K and Tc. The 
curved lines are a quadratic fit below 4 K. 




7^ (K^) 

FIG. 4: (color online). ps{T) plotted versus T^. The straight 
lines are quadratic fits to the data below 4 K except in the 
case of the lowest doping (Tc = 3 K), where the fit is to just 
below Tc. The data are linear in T^ up to T 5 K. 



approximation, giving a = ai— 1(72 = iujfio/Zg for the mi- 
crowave conductivity. The superfluid density is defined 
to be = 1/A^ = uoiiQ(j2' 

Broadband spectroscopy of the quasiparticle conduc- 
tivity (Ji(r^,T) has been carried out using bolometric 
measurements of RAVt^T) between 0.1 and 20 GHz, as 
described in Refs. ill and [H. The YBa2Cu306.333 el- 
lipsoid and a Ag:Au reference sample were positioned in 
symmetric locations at the end of a rectangular coaxial 
transmission line, with the microwave JT-field again ori- 
ented along the c axis of the ellipsoid. Rs{^^ T) has been 
inferred from the synchronous rise in sample tempera- 
ture in response to incident microwave fields modulated 
at 1 Hz. The Ag:Au sample acts a power meter, provid- 
ing an absolute calibration. At low frequencies, <Ji can 
be obtained from Rs from a knowledge of the penetra- 
tion depth: in this limit ai ~ 2Rs/Vfiii^\^. At higher 
frequencies, the quasiparticle conductivity starts to con- 
tribute to electromagnetic screening, effectively reduc- 
ing A. The quasiparticle shielding effect must be taken 
into account self-consistently, and the procedure for do- 
ing this is described in detail in Appendix C of Ref. [82|. 
As part of this process, the quasiparticle contribution to 
(72 is inferred from a Kramers-Kronig transform of cri (yt). 
This in turn requires a robust means of extrapolating 
cFi(yt) outside the measured frequency range. In previ- 
ous work^SiiS^ we have shown that the phenomenological 
form, 

am=a^l[l + {n/V)y], (6) 

works well for cuprate superconductors, with the expo- 
nent y ranging from 1.4 to 1.7. A Drude model, on the 
other hand, corresponds to ^ = 2. Physically, the non- 



Drude exponents stem from the strong energy depen- 
dence of scattering rate in an unconventional supercon- 
ductor. At low temperatures, thermally excited quasipar- 
ticles make a relatively small contribution to electromag- 
netic screening, so the extraction of (Ji(p) from Rs{^) 
is not particularly sensitive to variations in y. A similar 
procedure is used to estimate the quasiparticle conductiv- 
ity spectral weight: in that case there is more sensitivity 
to the choice of exponent when integrating (Ji{Q). 

V. RESULTS AND DISCUSSION 

psiT) is plotted in Fig. [3] for a subset of the dopings. 
The most prominent feature of the data is the linear 
T dependence of ps in the middle of the temperature 
range, which crosses over to a weaker temperature de- 
pendence at low T. The main questions about these 
data are: what is the limiting low temperature form of 
Ps{T)?; is the crossover the result of disorder?; and is 
the linear T dependence at higher temperatures charac- 
teristic of the behaviour of the ideal, clean system? To 
address these issues, we first look at the low temperature 
range in more detail. Fig. [4] plots the data from Fig. [3] 
vs. T^, showing that ps{T) indeed crosses over to accu- 
rately quadratic behaviour. For the lowest doping (the 
fully relaxed state with Tc 3 K), the sample has been 
remounted in our dilution refrigerator system and mea- 
sured down to T = 0.05 K. This data is plotted vs. T^ 
in Fig. [5l We see that the quadratic behaviour is robust 
to the lowest temperatures, neither flattening out to ac- 
tivated behaviour nor turning up to reveal a power law 
intermediate between and T^. 

To test whether the curvature is the result of disorder. 
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FIG. 5: (color online). For the lowest doping in this study 
(Tc = 3 K), Ps{T) has been measured down to T = 0.05 K. 
The data, plotted versus T^, reveal that the asymptotic low 
T behaviour is quadratic in temperature. 



we switch now to broadband microwave spectroscopy, 
which probes the spectral weight of the zero-energy quasi- 
particles. Fig. [6] shows Rs{^) at T = 1.7 K for the 
Tc = 15.6 K doping. This has been converted to conduc- 
tivity ai{Q) in Fig.[7J using the self-consistent procedure 
described in the previous section. As mentioned above, 
we use a phenomenological form to fit to the conductiv- 
ity: ai{n) = ao/[l + {^/T)y]. Spectra with y = 1.4 and 
y = 1.7 provide equally good fits to the Rs{^) data in 
Fig. [6] — a Drude fit {y = 2), however, shows marked 
deviations at the high frequency end. At low frequency 
there is a narrow peak in (Ji{Q)^ of uncertain origin, that 
may be a fluctuation effect. In any case we are content 
to omit it from the fitting procedure as it contains an in- 
significant fraction of the total oscillator strength. Using 
the phenomenological model of conductivity, we calculate 
the uncondensed spectral weight, for different choices of 
exponent. Expressed in superfluid density units, we ob- 
tain Aps = 1.05 /im~^ for y = 1.4 and Apg = 0.70 /im~^ 
for y = 1.7. Fig. [3] also shows linear and quadratic fits 
to Ps{T) at low temperature. For comparison with the 
integrated T = 1.7 K spectral weight in cri(l^), we should 
use the difference between the linear extrapolation of ps 
to T = 0, and ps{T = 1.7 K): this is Aps = 1.03 pm''^. 
As this falls within the range estimated from integrating 
cri(r^), we conclude that the crossover to behaviour in 
ps (T) is most likely a disorder effect in an otherwise pure 
dx2_y2 state, and that linear fits to Ps{T) in the middle 
of the temperature range should provide a good measure 
of the low temperature slope in the absence of disorder. 

A useful characterization of the strength of disorder is 
provided by the temperature at which ps{T) crosses 
over from quadratic to linear behaviour. Using an in- 
terpolation formula, Aps{T) = AT'^/{T + 2Trf), simi- 
lar to that of Ref. [tO, Td is defined to be the point 
at which the slope of the high temperature linear be- 
haviour, Aps = ceT, matches the slope of the low tem- 
perature quadratic behaviour, Aps = /3T^. Using values 




5 10 15 

Frequency (GHz) 



20 



FIG. 6: (color online). Broadband bolometric measurement 
of the surface resistance, i?s(Q), at T = 1.7 K. Data are for a 
doping state with Tc — 15.6 K. The solid line is a fit using the 
phenomenological conductivity model, Eq.[6l with y — 1.7. A 
fit with y — l.A IS practically indistinguishable and provides 
an equally good representation of the data. The dashed line, 
a best fit to the Drude model {y — 2)^ shows clear deviations 
at high frequencies. 
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FIG. 7: (color online). The real part of the conductivity 
spectrum determined from the Rs{p) data in Fig.[6l The solid 
line is a fit to the conductivity spectrum for y — 1.7 using the 
phenomenological model Eq. [6] The small, narrow peak at 
low frequencies is a robust result of the analysis and indicates 
long lived currents, possibly associated with superconducting 
fluctuations. 



of a and (3 obtained from fits similar to those shown 
in Figs. [3] and H we plot Td = a/2(3 in Fig. El The 
crossover temperature lies between 4 K and 5 K at these 
low dopings. This is larger than the cross-over tempera- 
ture in the best samples of Ortho-II YBa2Cu306.5o and 
Ortho-I YBa2Cu306.99, where Td is less than 1 K. This 
is consistent with the lower degree of CuO chain order 
in YBa2Cu3 06.333, which is known to be the dominant 
source of residual scattering in the best YBa2Cu306+2y 
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FIG. 8: (color online). The doping dependence of the dis- 
order crossover temperature Td and the uncondensed super- 
fluid density Aps. Td is the temperature at which linear 
and quadratic fits to ps{T) match in slope, as defined in the 
text. Ay9s is the uncondensed spectral weight predicted from 
the difference between linear and quadratic extrapolations of 
Ps{T) to T = 0, and is consistent with the residual conduc- 
tivity spectral weight directly measured at Tc — 15.6 K via 
broadband spectroscopy. 



of such effects. We show in Fig.[9]that the idxy state only 
becomes visible once A^^^^ > kBTd. Similarly, we would 
expect the clean-limit behaviour of the B// state to be 
apparent once 4Acc > ksTd^ meaning that if a perturba- 
tion of the form Eq. [4] is present, then Acc must be 1 K 
or less^^ The constraints become even tighter in Ortho- 
II YBa2Cu306.5o and Ortho-I YBa2Cu3 06.99, where the 
disorder scale Td is less than 1 K. 

Finally, while we can rule out nested spin and charge 
density waves, our data say very little about commensu- 
rate orders that connect parts of the Fermi surface away 
from the nodes, as these will generally not alter the low 
energy spectrum. One such a scenario has been revealed 
by recent STM measurements on Bi2Sr2CaCu208+(5r^ 
which show ordered, nondispersing modulations of the 
DOS at high energies and simultaneously, at low ener- 
gies, arcs of Bogoliubov quasiparticles associated with 
the nodal dx2-y2 spectrum. The 'Bogoliubov arcs' ap- 
pear to terminate on the Bragg plane joining (0, tt) and 
(7r,0) points, leaving the nodal spectrum intact. This 
would be compatible with the conclusions we draw here 
about YBa2Cu306+2;. 



VI. CONCLUSIONS 



samples^ and is likely enhanced by proximity to the 
Mott insulator. Also shown in Fig.[8]is the residual DOS, 
expressed in superfluid density units as Aps, and inferred 
from the difference between linear and quadratic extrap- 
olations of Ps{T) to T = 0. Aps falls on underdoping, 
but remains a roughly constant fraction of ps{T = 0), 
consistent with the weak doping dependence of T^. 

We are able to draw tight conclusions from these mea- 
surements about the types and magnitudes of electronic 
order than might be competing with pure dx2-y2 su- 
perconductivity in YBa2Cu306+2y. We emphasize that 
to do this it is essential to have measurements of both 
the asymptotic low temperature form of Ps(T), and the 
residual DOS from ai{Q). On the basis of the limiting 
quadratic T dependence, which we have followed down to 
0.05 K, we can rule out any of the clean-limit behaviours 
shown in Table [U as well as the dx2_y2 + is state in the 
presence of disorder. We can also exclude the BCS-BEC 
crossover scenario, which predicts a T^/^ term in Ps{T) 
from incoherent Cooper pairs excited from the conden- 
sate. When disorder is included, four of the remaining 
states in Table |T] are compatible with quadratic behaviour 
in Ps{T). Of these, nested spin and charge density waves 
can immediately be eliminated, as they are not expected 
to be accompanied by a residual DOS. Of the remaining 
three, the simplest possibility is pure dj.2_y2 supercon- 
ductivity in the presence of a small amount of strong 
scattering disorder. However, we cannot rule out a small 
idxy component, nor a weak 0//-type circulating current 
phase. Nevertheless, we can place tight limits on the size 



We have shown that measurements of superfluid den- 
sity can be used as a sensitive probe of electronic orders 
than might compete with pure dx2-y2 superconductivity. 
Broadband conductivity measurements provide comple- 
mentary information on zero-energy quasiparticles that 
would be difficult to infer from ps{T) alone. Measure- 
ments on underdoped YBa2Cu306.333 reveal a crossover 
from linear to quadratic behaviour in ps (T) below a tem- 
perature 4 K to 5 K. The power law has been 
followed as low as 0.05 K and appears to be the asymp- 
totic low temperature behaviour. It is also accompanied 
by a residual quasiparticle spectral weight of correspond- 
ing magnitude, leading us to conclude that the crossover 
is a disorder effect. The observations immediately allow 
us to rule out BCS-BEC crossover physics; competition 
from dx2-y2 + is superconductivity; and spin and charge 
density waves that nest the nodal points. Due to the 
presence of disorder, we cannot eliminate the possibil- 
ity of either disordered dx2-y2 + idxy superconductivity, 
provided Ad^y < 4 - 5 K; or a perturbation of the form 
Eq. [Hfrom a 6//-type circulating current phase, as long 
as Acc ^ 1 K. The small magnitude of the term is com- 
patible with related observations from /iSR^^ neutron 
scattering^ and polar Kerr-effect measurements 
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APPENDIX A: d^2_y2, d^2_y2 +idxy AND d^2_y2 + is 
STATES 



The dj.2_y2 



idxy state and the dj.2_y2 



\s state are 

two candidate order parameters that may compete with 
pure dx2_y2 superconductivity in the cuprates. In this 
appendix we review the theory of the penetration depth 
in the presence of disorder and gauge the extent to which 
these states can be distinguished by microwave experi- 
ments. The theory of unconventional superconductivity 
in the presence of elastic scattering disorder has been de- 
veloped by many authors ^i^i^^SiiSS^^^^ j^^g 

been reviewed in several places.— i^*^ In these systems, 
disorder not only imparts a finite lifetime to the quasipar- 
ticles, it alters the excitation spectrum by pair-breaking, 
and the two effects must be dealt with together. The self- 
consistent t-matrix approximation (SCTMA) provides a 
powerful approach for capturing this physics, particularly 
in the resonant scattering limit, where the impurity is on 
the verge of binding a quasiparticle at the Fermi energy. 
In the SCTMA, impurities are usually approximated as 
point defects that scatter in the 5-wave channel. The 
effect of the disorder is to renormalize the quasiparticle 
energy uo and the superconducting gap Ak, which can be 
expressed in the following way: 



uo ^ Co = uo ^ \7lT- 



N{uj) 



Ak + iTrr 



' + Ar2(cj)+P2(cj) 



(Al) 



. (A2) 



Here T = nin/7r'^D{eF), where rii is the impurity concen- 
tration, n is the conduction electron density, and D{eF) 
is the density of states at the Fermi level. ^- The impurity 
scattering strength is characterized by c, the cotangent 
of the s-wave scattering phase shift. The quasiparticle 
density N{uo) and pair density P{uo) depend on details of 
the particular superconducting state and are defined be- 
low for the different types of order parameter. For purely 
unconventional order parameters, (Ak)FS = and P{uo) 
vanishes — these states are therefore unrenormalized by 
s-wave scatterers. 

We are primarily interested in the behaviour of the low 
energy excitations so, without loss of generality, we take 
the two-dimensional Fermi surface to be isotropic, and 
the gap functions to be the simplest cylindrical harmonics 
of the required symmetry: 



^d^2_y2 = Aocos20 , 
Ad,^ = 7/Aosin2(/) . 
A, = CAo . 



(A3) 
(A4) 
(A5) 



Here (j) measures angle from the Cu-0 bond direction and 
T] and C are constants. For the pure dj.2_y2 state there is 



no gap renormalization. The quasiparticle density is 
N{u;) = 



uo 



= -k(^) , (A6) 



y^Uj'^ - Al C0S2 2(1)/ 7T 



where (...)(/, is an angle average around the cylindrical 
Fermi surface, K{x) is the complete elliptic integral of 
the first kind, and the branch of the square root in Eq. lAll 
is chosen so that uo has positive imaginary part. In the 
strong-scattering (unitarity) limit, for instance, c = 
and uj{uj) is a root of 

u = uj^ 'm^r/2K{Al/Cj'^). (A7) 

uo{uo) encodes all the physics of scattering and pair- 
breaking. Inserted into the real part of Eq. lA6l it gives the 
quasiparticle density of states in the presence of disorder. 
To calculate penetration depth using a), a modification of 
Eq.[2]is used>^ 



2fcBT \(^2_A2)i 



k'' / FS 



The density of states factor is 



Al cos2 : 



FS 



(a;2- Ag cos2 2(^)i 



r,2 



(A8) 



(A9) 



Al-oj^^ 



where E{x) is the complete elliptic integral of the second 
kind. 

For the dx2_y2 + idxy state, A((/)) = Ao(cos2(/) + 
i?7sin2(/)), and there is similarly no gap renormalization. 
The quasiparticle density is 



N{uj) 




uo 



Ag(cOs2 20 + 7^2 gij^2 20)^ 

'(l-^')Ar 



(AlO) 



TT ^CO^ - 7^2 A2 



K 



The density of states factor in Eq. IA8I becomes 



A^ 



(uo^ - Al)i / 

I ^ k W FS 



/ Ag(cos2 2(/) + 7/2sin^2(/)) 
\ (a;2 - A2(cos2 2^ + 7y2 sin^ 20)) ■ 
2 1 

7^ ^UO^ - 7^2 a| ^ 

(l-^^)Ag 



(All) 



K 



a;2 - 7^2 A2 



A2 
^0 









-r,^Al)\ 
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FIG. 9: (color online). The onset of quadratic temperature 
dependence of ps{T) in a (1^2 _y2 + idxy superconductor as 
a function of disorder strength, relative to that of a (1^2 _y2 
state. [3d+id is the coefficient of ps(T) for the (1^2 _y2 -\-idxy 
superconductor. jSa is the same quantity for the d^2_y2 state. 
Disorder strength is characterized by the disorder crossover 
temperature Td of the d^2_y2 superconductor, as defined in 
the text. Data are plotted for different values of the dxy gap, 
Ad^^, and scale well as a function of Td/ Ad^y. The two pair- 
ing states are difficult to distinguish on the basis of Aps{T) 
once ksTd > A^^^. 



In the dj.2_y2 + \s state, impurity renormalization of Ag 
must be taken into account. The renormalization equa- 
tions lAll and IA2I can be rewritten 



1 = ^ + 

A,, 



iTrr 



N{lu)/Co 



iTrr 



C2 + A/-2(w) +P2(^) 



where 




+N^{oj)+P^uj) 



Ag cos2 20 - A2 , 



(D2 - A2 cos2 20 - A2 , 



(A12) 
(A13) 

(A14) 
(A15) 



= P{(jo)/Ks^ the quantities uo/Cb and 
A5/A5 obey identical equations and therefore A5 = 
AsCo/uo. IA12l and lA13l can then be combined into a single 
equation 



a; = cj + IttT 



c2+A/-2(cj)(l+A2/cj2) ' 



(A16) 



where 



N{uj) 



\ ^uj^il - /S?Juj^) - Ag cos2 20 ^ 
2 1 



A^ 
^0 



The corresponding term in Eq. 



(A17) 



IS 



(w2 



AD* 



FS 



Ag(cos2 2(j)^rf sin^ 20) 



A2(cos2 20- 

2 1 



7^2 sin^ 20))- 



TT ~ 



i_ 

J- , .9 



K\ 



A2 
^0 



(^2 



.(1- 



A2n 



Ai-(1- 



A2^ 



-E\ 



^0 



.(1- 



(A18) 



We are now in a position to compare results for the 
three order parameters. The forms for the density of 
states Niuo) and the superfluid density ps{T) are shown 
in Table [U both in the clean limit and in the presence 
of strong scattering disorder (c = 0). The key feature of 
the clean dx2_y2 + idxy and d^-^-y^ + \s states is a finite 
energy gap, giving rise to activated behaviour in ps{T). 

:2_y2 -\-is states behave very dif- 



The d^2_ 



- idxy and d, 



ferently in response to disorder. In the dx2-y2 + is case, 
the energy gap is robust. This can be traced back to the 
expressions for the renormalized frequency, Eqs. lA16l and 
IA171 Impurity renormalization of A5 leads to solutions 
for LJ that are purely real for cj < A5, preventing the for- 
mation of any low- lying quasiparticle states in N{uj)^ 
The dx^-y^ + idxy case is quite different: pair breaking 
occurs for even small amounts of disorder, leading im- 
mediately to a term in ps{T). The term starts 
out weak, but grows in magnitude until it is comparable 
to that of a pure dx2-y2 superconductor with a similar 
amount of disorder. This cross over is charted in Fig. [9l 
which shows that the dx^-y^ and dx2-y2 + idxy states 
become indistinguishable when the energy scale for the 
disorder, ksTd^ becomes comparable to A^^^ . 
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